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ABSTRACT 

We have discovered a bow shock shaped mid-infrared excess region in front of 5 Velorum using 
24 /jm observations obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The excess 
has been classified as a debris disk from previous infrared observations. Although the bow shock 
morphology was only detected in the 24 /urn observations, its excess was also resolved at 70 /im. 
We show that the stellar heating of an ambient interstellar medium (ISM) cloud can produce the 
measured flux. Since S Velorum was classified as a debris disk star previously, our discovery may call 
into question the same classification of other stars. We model the interaction of the star and ISM, 
producing images that show the same geometry and surface brightness as is observed. The modeled 
ISM is ~ 15 times overdense relative to the average Local Bubble value, which is surprising considering 
the close proximity (24 pc) of 5 Velorum. 

The abundance anomalies of A Bootis stars have been previously explained as arising from the same 
type of interaction of stars with the ISM. Low resolution optical spectra of 6 Velorum show that it 
does not belong to this stellar class. The star therefore is an interesting testbed for the ISM accretion 
theory of the A Bootis phenomenon. 

Subject headings: stars: evolution - stars: imaging - stars: individual (HD 74956, 6 Velorum) - 
ISM: kinematics and dynamics - infrared: ISM - radiation mechanisms: thermal, 
Shockwaves 



1. INTRODUCTION 

Using IRAS data, more than a hundred main-sequence 
stars have been found to have excess emission in the 12 - 
100 //m spectral range (jBeckman fc Parescel 19931 ). Many 
additional examples have been discovered with ISO and 
Spitzer. In most cases the spectral energy distributions 
(SEDs) can be fitted by models of circumstellar debris 
systems of thermally radiating dust grains with temper- 
atures of 50 to 200 K. Such grains have short lifetimes 
around stars: they either get ground down into tiny dust 
particles that are then ejected by radiation pressure, or if 
their number density is low they are brought into the star 
by Poynting- Robertson drag. Since excesses are observed 
around stars that are much older than the time scale for 
these clearing mechanisms, it is necessary that the dust 
be replenished through collisions between planetesimals 
and the resulting collisional cascades of the products of 
these events both with themselves and with other bod- 
ies. Thus, planetary debris disks are a means to study 
processes occurring in hundreds of neighboring planetary 
systems. Spitzer observations are revealing a general re- 
semblance in evolutionary time scales and other proper- 
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ties to the events hypothesized to have occurred in the 
early Solar System. 

Although the planetary debris disk hypothesis appears 
to account for a large majority of the far infrared ex- 
cesses around main-sequence stars, there are two alter- 
native possibilities. The first is that very hot gas a round 
the stars is res ponsi ble for free- free emission (e.g., ICotd 
U987tlSu et"aTH2006fh The second possibility is that the 
excesses arise through heating of dust grains in the in- 
terstellar medium around the st ar, but not i n a bo und 
structure such as a debris disk. I K alas et al.l (|2002l ) no- 
ticed optical reflection nebulosities around a number of 
stars with Vega-like excesses. These nebulosities show 
asymmetries that would not be typical of disks, they 
have complex, often striated structures that are remi- 
niscent of the Pleiades reflection nebulosities, and they 
are much too lar ge in extent to be gr avitationally bound 
to the stars fsee lGorlova et a l. 2006). 

Dynamical rather than st ationary interactions with the 
ISM are more interesting (jCharbonneaul fl991h . Origi- 
nally, it was proposed that ISM dust grains could interact 
direc t ly with material in deb ris disks ( Lis sauer fc G riffith 
119891 : IWhitmire et al.l fl992). However, it was soon re- 
alized that photon pressure from the star would re- 
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TABLE 1 

The parameters of <5 Velorum 



F 2i * 


^70* 


PISM 




-Fstar24^ 




^excess24^ 


-?cxccss70 


(mjy) 


(mjy) 


(10" 24 g cm- 3 ) 


(km s _1 ) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


1569 ± 42 


237 ± 50 


5.8 ±0.4 


36 ±4 


1277 


147 


174 


141 



* Observed fluxes with the large aperture 

T Photospheric values - not including G star component 

t Modeled excesses at large aperture 



pel interstellar grains, resulting in grain-free zones with 
possible bow-shock geometry around luminous stars 
(|Artvmowicz fc Clampidll9971h 

This scenario has been proposed to account for the 
abundance anomalies associated with A Bootis stars. 
These are late B to early F-type, Population I stars with 
surface underabundances of most Fe-peak elements and 
solar abundances of lighter elements, su ch as C, N, O and 
S. In the diffusion /accretion model (|Venn fc Lambertl 
Il990t iKamp fc Paunzenl 120021: iPaunzen et al.ll2003f) . it 
is suggested that the abundance anomaly occurs when 
a star passes through a diffuse interstellar cloud. The 
radiation pressure repels the grains, and hence much of 
the general ISM metals, while the gas is accreted onto 
the stellar surface. While the star is within the cloud, a 
mid-infrared excess will result from the heating of the in- 
terstellar dust; however, after the star has left the cloud 
the abundance anomalies may persist for ~ 10 6 yr in its 
surface layers (|Turcotte fc Charbonneaulll993D without 
an accompanying infrared excess. 

There have been few opportunities to test the predic- 
tions for dynamical interactions of main-s equence stars 
with the ambient interstellar medium. iFrance et al.l 
(|2007f ) have studied a bo w shock generated by the 09.5 
runaway star HD 34078. lUeta et alj (|2006l ) describe the 
bow shock betwee n the mass loss wind of the A GB star R 
Hya and the ISM. iNoriega-Crespo et all (|1997ft identified 
58 runaway OB stars with an observable bow shock struc- 
ture using high res olution IRAS 60 /im emission maps. 
iRebull et alj (|2007f) discovered that the young B5 star 
HD 281159 is interacting with the ISM, producing spher- 
ical shells of extended IR emission centered on the star 
with a spike feature pointing from the star into the shells. 
None of these cases correspond to the type of situation 
that might be mistaken for a debris disk, nor which would 
be expected to produce a A Bootis abundance pattern. 

S Velorum is a nearby (~ 24 pc) stellar system (at least 
five members) 1 , with modest excess in the IRAS data. It 
has been cla ssified as an A- type star with a debris disk 
syste m (e.B:.JAumajid [l 9^[l9^[Cotdll987t lC^ii et all 
I2006t ISu et all l2006fh lOtero et al.l (|2000l ) observed a 
drop in the primary component's brightness (~ 0™3) and 
showed that it is an eclipsing binary with probably two 
A spectral type co mponents. With the available data, 
lArgyle et al.l (|2002) computed the system's parameters. 
They suggested that the eclipsing binary (Aa) consists 
of two A dwarfs with spectral types A1V and A5V and 
masses of 2.7 and 2.0 M Q and with separation of 10 mas. 
The nearby B component is a G dwarf with mass around 



It is a c omple x multiple system: I Otero et alj ( 200C ) 
Hanburv et~ap JT974T); IHorch et alj ||2000T ); lArgyle et all (120021 ) 
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1 M and separation of 0'.'6 from the main component. 
There is also another binary (CD component) at 78" from 
the star. 

In §2, we report measurements demonstrating that 
this star is producing a bow shock as it moves 

through an interstella r cloud as hypothesized by 

lArtvmowicz fc Clampinl (|1997h . In §3, we model this be- 
havior using simple dust grain parameters and show sat- 
isfactory agreement with expectations for the ISM and 
properties of the star. We discuss these results in §4, 
where we show that the star is most likely not part of 
the A Bootis stellar class. Thus, 5 Velorum provides a 
test of the diffusion/ accretion hypothesis for A Bootis 
behavior. 

2. OBSERVATIONS AND DATA REDUCTION 

We present observations of 6 Velorum at 24 and 70 fim 
obtained with the Multiband Imaging Photometer for 
Spitzer (MIPS) as part of three programs: PID 57 (2004 
Feb 21), PID 20296 (2006 Feb 22, Apr 3) and PID 30566 
(2006 June 12). For PID 57, we used 3 second expo- 
sures at four dither positions with a total integration 
time of 193 seconds. The other observations at 24 /im 
(PID 20296) were done in standard photometry mode 
with 4 cycles at 5 sub-pixel-offset cluster positions and 3 
sec integrations, resulting in a total integration of 902 sec 
on source for each of the two epochs. The star HD 217382 
was observed as a PSF standard (AOR ID 6627584) for 
PID 57, with the same observational parameters. The 
observation at 70 /im (PID 30566) was done in standard 
photometry default-scale mode with 10 sec integrations 
and 3 cycles, resulting in a total integration of 335 sec 
on source. 

The binary component Aa was not in ecl ipsing phase 
accord ing to the ephemeris equations by lOtero et al.l 
(2000) at either epoch. The period of the eclipse is 
~ 45.16 days, and the system was ~ 13 days before a 
primary minimum at the first, ~ 3 days before one at 
the second and ~ 7.7 days before one at the third epoch 
for the 24 /im observations. The 70 /zm observation was 
2.53 days before a secondary minimum. 

The data were processed usin g the MIPS instrum ent 
team Data Analysis Tool (DAT. iGordon et al.ll2005f) as 
descri bed bv lEngelbracht et alj (|2007l ) and IGordon et al.l 
(p00l . Care was taken to minimize instrumental arti- 
facts (details will be discussed in an upcoming paper, Su 
et al. 2007, in preparation). 

Fitting the model described later demands flux mea- 
surements within a constant large external radius (see 
details in §3). Therefore, photometry for the target was 
extracted using aperture photometry with a single aper- 
ture setting. The center for the aperture photometry at 
both 24 and 70 ^m was determined by fitting and cen- 
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Fig. 1. — The panels show 24 /im images of S Velorum. All images are in logarithmic scaling, the FOV is ~ 2.' 74 X 2'34. The scaling 
of the images are: —0.5 — 4 MJy sr — 1 . Top-left panel: The original observed composite image from the 2 nd and 3 rd epochs. Top-right 
panel: PSF oversubtracted image, which shows the bow shock structure far from the star. Bottom-left panel: The intensity scaled PSF 
subtracted image (first epoch), which shows the bow shock structure close to the star. This image shows the orientation of the images and 
the proper motion direction of the star. The arrow bisecting the bow shock contour shows the calculated direction of the modeled relative 
velocity. Bottom-right panel: Same image as the bottom-left panel, but with intensity contours plotted. The intensity contours are at 0.25, 
1.0, 1.75, 2.5 and 3.25 MJy sr — 1 from the faintest to the brightest, respectively. The contours show that the extended emission consists of 
incomplete spherical shells, centered on S Velorum. 

troiding a 2-D Gaussian core. A radius of 56'.'025 was 
used for both wavelengths, with sky annulus between 
68'.'95 and 76'.'34. The aperture size was chosen to be 
large enough to contain most of the flux from the bow 
shock, but small enough to exclude the CD component 
to avoid contamination. The CD component was bright 
at 24 /im at a distance of 78" from the AaB compo- 
nents, but could not be detected at 70 /im. Aperture 
corrections were not applied because of the large size of 
the aperture. Conversion factors of 1.068 x 10 -3 and 
1.652 x 10 1 mJy arcsec^ 2 MIPS.UNIT -1 were used to 
transfer measured instrumental units to physical units 
at 24 and 70 /im, respectively. 

Faint extended asymmetric nebulosity offset from the 
central star is apparent at 24 /im, with the dark Airy 
rings partially filled in. Using standard aperture and 
point-spread-function (PSF) fitting photometry opti- 
mized for a point source, the total flux is 1420±42 mJy, ~ 
1.12 times the expected photospheric flux, which was de- 



termined by fitting a Kurucz model (|Castelli fc Kuruc3 
l2003f ) to the optical and near infrared photometry and 
extrapolating it to 24 and 70 /im. The large aperture 
photometry value is greater by another factor of ~ 1.1, 
which puts it above the expected photospheric flux by 
a factor of ~ 1.25. The final photometry measurements 
(using the large aperture setting) are listed in Table [TJ 
We also list the modeled photospheric flux of the star and 
the modeled value of the IR excess. Since the measured 
excess depends on the aperture used, to avoid confusion 
we do not give a measured excess value, only the pho- 
tospheric flux which can be subtracted from any later 
measurements. The photospheric flux given in Table [T] 
does not include the contribution from the G dwarf (90 
and 10 mJy at 24 and 70 /mi, respectively). The top left 
panel in Figure 1 shows the summed image from epochs 
2 and 3, to demonstrate the asymmetry suggested even 
before PSF subtraction. 
For the first epoch 24 /im image, the reference star im- 
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age was subtracted from the image of S Velorum, with 
a scale factor chosen as the maximum value that would 
completely remove the image core without creating sig- 
nificant negative flux residuals. The deeper exposures 
from the second and third epochs were designed to re- 
veal faint structures far from the star, where the observed 
PSF is difficult to extract accurate ly. Therefor e, we used 
simulated PSFs (from STinyTim (Kristll2002D ) and the 
MIPS simulator 2 . Because bright structures nearly in 
the PSF contribute to the residuals at large distances, we 
oversubtracted the PSF to compensate. The first epoch 
PSF subtracted 24 fim image is shown in the bottom 
panels of Figure [T] and the composite from epochs 2 and 
3 in the upper right. 

The PSF subtracted images in Figure [T] show that the 
asymmetry is caused by a bow shock. As shown in the 
lower left, the head of the bow shock points approxi- 
mately toward the direction of the stellar proper motion. 
The bottom right panel shows the excess flux contours 
and that it consists of incomplete spherical shells cen- 
tered on 6 Velorum. Combined with the upper right 
image, there is also a parabolic cavity, as expected for 
a bow shock. The stagnation points (where photon pres- 
sure equals gravitational force) of the grains in the bow 
shock are within ~ 200 AU of the star, according to 
the observations. A notable feature in the upper right 
is the wings of the bow shock, which are detectable to 
- 1500 AU. 




Fig. 2. — The panels show the 70 (im image of 5 Velorum. All 
images are scaled logarithmically from —0.5 - 3 MJy sr _1 . The 
FOV is ~ 2.' 46 X 3'03. The orientation of the images is the same 
as in Figure [T] First panel: the observed image. Middle panel: 
the PSF subtracted image. The residual flux seems close to being 
concentric. Last panel: the intensity contours. They suggest that 
there is a faint concentric 70 fim excess further from the star that 
fades at the cavity region behind the star. 



The 70 fim observation is shown in Figure [2l The PSF 
subtraction (scaled to the point source flux of 125 mJy) 
does not reveal the bow shock structure at this wave- 
length, only that there is extended excess. The to- 
tal flux of the residual of the PSF subtracted image is 
119 mJy. The intensity contours (last panel) suggest that 
the 70 /jm excess fades at the cavity behind the star, but 
the effect is small. The geometry and direction of the 
bow shock are discussed in more detail in §3.2. 

3. THE BOW SHOCK MODEL 

Bas ed on a previous su g gestio n by IVenn fc Lambertl 
(1990), iKamp fc Paunzer] (|2002f ) proposed a physical 
model to explain the abundance pattern of A Bootis stars 

2 Software designed to simulate MIPS data, including optical 
distortions, using the same observing templates used in flight. 



through star-ISM interaction and the diffusion/accretion 
hypothesis. Their model is based on a luminous main- 
sequence star passing through a diffuse ISM cloud. The 
star blows the interstellar dust grains away by its ra- 
diation pressure, but accretes the interstellar gas onto 
its surface, thus establishing a thin surface layer with 
abundance anomalies. So long as the star is inside the 
cloud, the dust grains are heated to produce excess in 
the infrared above the photospheric radiation of the star. 
Martmez-Galarza et al. (2007, in prep.) have developed 
a model of this process and show that the global spectral 
energy distributions of a group of A Bootis type stars that 
have infrared excesses are consistent with the emission 
from the hypothesized ISM cloud. Details of the model 
can be found in their paper. Here we adapt their model 
and improve its fidelity (e.g., with higher resolution in- 
tegrations), and also model the surface brightness distri- 
bution to describe the observed bow shock seen around 
S Velorum. 

3.1. Physical description of the model 

The phenomenon of star-IS M interactions generating 
bow s hocks was first studied by I Artvmowicz fc Clampinl 
(|1997| ). They showed that the radiative pressure force 
on a sub-micron dust grain can be many times that of 
the gravitational force as it approaches the star. The 
scattering surface will be a parabola with the star at the 
focus point of the parabolic shaped dust cavity. Since 
the star heats the grains outside of the cavity and close 
to the parabolic surface, an infrared-emitting bow shock 
feature is expected. 

The shape of the parabola (for each grain size) 
can be given in terms of the distance between the 
star (focus) and the vertex. This so-called avoid- 
ance radius (or the p/2 parameter of the scattering 
par abola) can be calculated from energy conservation to 
be (jArtymowicz fc Clampinlll997f ): 

_ 2(/3°-l)GM 

'av 9 ' V-V 

where a is the radius of the particle, M is the mass of 
the star and v re i is the relative velocity between the star 
and the dust grains. 

P a is the ratio of photon pre ssure to gravitatio nal force 
on a grain and it is given by ([Burns et al.lll979l ): 

P a = 0-57Qp t -Uwr (— y 1 ( -f-s) " . ( 2 ) 

^ pr M/M© \nmj Vgcm y 
where 6 is the bulk density of the grain material and Q® T 
is the radiation pressure efficiency averaged over the stel- 
lar spectrum. Qp r (A) can be expressed in terms of grain 
properties (absorption coefficient Q% b (\), scattering co- 
efficient Q^ na (X ) and the scatter i ng asymmetry factor 
q = (cos a (A)), (Burns et alj|197ft iHenvev fc Greensteinl 
Il93l : 

Qpr(A) = Q: b (A) + Q: ca (l- 5 ), (3) 

which gives 

_ / Qg r (A)-B(r„, A)dA 
^ pr fB(n,X)dX ' [) 

where B(T*, A) is the Planck function. We adopted astro- 
nomical silicates in our model with 8 = 3.3 g cm -3 from 
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Fig. 4. — The panels show the 24 fim morphology of the bow shock viewed at different inclinations, starting from 90° (left), 70° (middle) 
and 50° (right). 



iDraine fe Lee] (119841) andlLaor fe Drairiel (Il993f). W e con- 
sidered a MRN (jMathis. Rumpl fe Nordsiecklll977ft grain 
size distribution in our model: 



dn = Ca 7 da, 



(5) 



where C is a scaling constant and n is the number density 
of the cloud with 7 = 3.5 and grain sizes ranging from 
0.005 - 0.25 ^. 

With these equations we are able to model the avoid- 
ance cavity for a grain that encounters a star with known 
mass, luminosity and relative velocity. The model de- 
scribes a situation where the expelled grains are instantly 
removed from the system rather than drifting away, but 
this only causes a minor discrepancy in the wing and 
almost none in the apex of the parabola compared to 
the actual case. In the actual scenario only those parti- 
cles get scattered back upstream that encounter the cen- 
tral star with small impact parameter (~ r" v /2). This 
means that most of the grains will get expelled toward 
the wings, where the grains go further out and emit less 
infrared excess, thus their contribution to the total flux 
will be small. 

The model determines the number density of certain 
grain sizes and the position of their parabolic avoidance 
cavity. Outside of the cavity we assumed a constant num- 
ber density distribution for each grain size. To calculate 
the surface brightness of the system and its SED we as- 
sumed a thermal equilibrium condition, with wavelength 
dependent absorption and an optically thin cloud. 

3.2. Model Geometry and Parameters 

The model described in §3.1 gives the distribution and 
temperature for each grain size. This model was imple- 
mented in two ANSI C programs. The first program 
fits the SED of the system to the observed photometry 
points, while the second program calculates the surface 
brightness of the system. The fitted photometry included 
uvby, UBV, HIPPARCOS V band, 2MASS, IRAS and 
MIPS (24 and 70 fim) data. We subtracted the 24 and 
70 fim flux contributed by the G star (90 and 10 mJy, 
respectively) from the MIPS observations, because we 
wanted to model the system consisting of the two A stars 
and the bow shock. 

The input parameters are: stellar radius, mass-to- 
luminosity ratio (MLR), relative velocity of cloud and 
star, ISM dust density, cloud external radius and the 




Fig. 3. — The nomenclature of the angles of the system. The 
heavy line is the grain avoidance parabola, if is the rotation angle 
of the system on the plane of the sky (our initial guess was 4° N 
from the calculated direction of relative motion shown in Figured}, 
i is the inclination and r av is an avoidance radius. The observer is 
viewing from the axis pointing to the bottom left. 



distance of the system. The stellar radius, MLR and the 
distance can be constr ained easily. We determi ned the 
best-fit Kurucz model (jCastelli fc Kuruc3 12003[ ) by fit- 
ting the photometry points at wavelengths shorter than 
10 fim. Since the distance is known to high accuracy from 
HIPPARCOS we can determine the radius and thus the 
lumin osity of the star pr ecisely. The mass was adopted 
from lArgvle et all (|2002D . The G dwarf's luminosity is 
only 1% of the system, so leaving the star out does not 
cause any inconsistency. Its mass is only 17% of the 
total mass, which can only cause minor changes in the 
determined final relative velocity, but none in the final 
surface brightness or the computed ISM density. The 
model then has three variable parameters: the density 
of the ISM grains (pdust - does not include gas), the rel- 
ative velocity between the cloud and the star (u re i) and 
the external radius of the cloud (r ext ). The model should 
describe the total flux from exactly the area used for our 
photometry. The aperture radius of 56'.'025 (1366 AU at 
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Fig. 5. — Left panel: The \ 2 phase space for Pdust vs - 
(left panel) cut at u rel = 36 km s . 



with constrained r cx t = 1366 AU. Right panel: The p vs. v re \ phase space 



the stellar distance of 24.45 pc) was used as r oxt . Both 
programs calculate the Qp r , f3 a , r° Y , n values and then 
the temperature at r° v for each grain size. 

The SED modeling program decreases the temperature 
value from the one at r av by 0.01 K steps and finds the 
radius for that corresponding grain temperature. The 
program does not include geometrical parameters such 
as the inclination or the rotation angle of the system, 
since these are irrelevant in calculating the total flux. It 
calculates the contribution to the emitting flux for every 
grain size from every shell to an external radius (r oxt ) 
and adds them up according to wavelength. 

The program that calculates the surface brightness 
uses a similar algorithm as the SED program, but it cal- 
culates the temperature at 1 AU distance steps from r" v 
for every grain size and calculates the total flux in the 
line of sight in I AU 2 resolution elements. 

The total inclination t of the bow shock was not in- 
cluded as a parameter, since by eye the observed images 
seemed to show an inclination of l k, 90° (a schematic 
plot of the angle nomenclature is shown in Figure [3|) . 
This approximation is strengthened by the radial veloc- 
ity of the star, which is only ~ 2 km s _1 compared to 
the tangential velocity of ~ 13 km s~ . This assures that 
the motion of the system is close to perpendicular to the 
line of sight. However, we have found that the bow shock 
has similar appearance for a significant range of angles 
(±20°) relative to l = 90°. We illustrate this in Figured 
If the relative velocity vector would have a 70° (or 110°) 
inclination it would only cause minor differences in the 
modeled velocity (Aw rc i w 3 km s' 1 ) and ISM density 
(ApisM « 0.2 x 10~ 24 g cm' 3 ). At an inclination of 50°, 
the "wings" spread out and the bright rim at the apex 
starts to become thin. 

With interst e llar Fell and Mgll measurements 
iLallement et al.l (|1995f ) showed that the Local Interstel- 



lar Cloud (LIC) has a heliocentric velocity of 26 km s 1 
moving towards the galactic coordinate /// = 186 ± 3°, 
b u = -16 ± 3°. Since S Velorum is at l n « 272°, 
bu Rs — 7°, the LIC is also moving perpendicular to our 
line of sight at the star and in the direction needed to 
reach a high rel a tive ve locity between the star and cloud. 
iCrawford et alj (|1998f ) showed a low velocity interstel- 
lar Ca K line component in the star's spectrum with 
^hoiio = 1-3 ± 0.4 km s _1 , which also proves that the 
ISM's motion is perpendicular to our line of sight at S 
Velorum. The offset of the proper motion direction of 
the star from the head direction of the bow shock by a 
few degrees could be explained by the ISM velocity. A 
simple vectorial summation of the star and the ISM ve- 
locities should give a net motion in the direction of the 
bow shock. 

3.3. Results 

We first tried to find the best fitting SED to the pho- 
tometry points corresponding to wavelengths larger than 
10 /Ltm (MIPS, IRAS) with \ 2 minimization in the w rc i 
vs. pdust phase space. We defined x 2 & s: 



X 



(Fo 



bs 



calc ) 



(6) 



J obs 



The x 2 phase space with r cx t = 1366 AU showed no 
minimum (Figure [5l left panel) . The interpretation of 
the diagram is as follows: if the relative velocity is small, 
then the avoidance radius will be large. Consequently 
the grains will be at relatively low temperature and the 
amount of dust required to produce the observed flux 
increases. On the other hand, if the relative velocity is 
large, then the grains can approach closer to the star and 
heat up to higher temperatures. As a result a smaller 
dust density is enough to produce the observed flux. 
Therefore, the combination of the density of the cloud 
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Fig. 7. — Left panel: \ m the phase space of ip vs. u re i. Right panel: The phase space cut at tp = —4°, showing the best fit for u re i. 
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Fig. 6.— The value of function of grain size. The solid 

lines are curves for a silicate bulk density of 3.3 g cm -3 , while the 
dashed ones are for 2.2 g cm -3 . The curves are for v Ic i values of 
25, 30, 35, 40 and 45 km s — 1 from top to bottom, respectively. 



and the relative velocity can be well constrained by the 
broad-band SED alone, but not each separately. 

By using surface brightness values from the observa- 
tions and the model calculations we were able to deter- 
mine the v Te \ parameter and thus eliminate the degener- 
acy of the model. Since the bow shock is a parabolic fea- 
ture it has only one variable, the avoidance radius (r^ ¥ ), 
which is the same as the p/2 parameter of the parabola 
(with p being the distance between the focus point and 
the vertex). The value of r° v does change as a function 
of grain size, but the head of the bow shock will be near 
the value where the avoidance radius has its maximum 
as a function of grain size. As can be seen in Figure 
the avoidance radius has a maximum at ~ 0.06 /im grain 
size. The value of the avoidance radius on the other 



hand only depends on the relative velocity between the 
ISM cloud and the star. This way we can constrain the 
second parameter of the model (v TC \). The relative ve- 
locity has to be set so that the avoidance radius of the 
~ 0.06 /im grain is around half the parabola parameter 
value. This method gives a value that only approximates 
the true one, but it can be used as an initial guess. 

The i> ro i parameter was constrained by comparing 
the PSF subtracted image "wings" with model im- 
ages. Within a range of ±6 km s _1 of our initial guess 
(vrd = 35 km s _1 ) with 1 km s^ 1 steps, we generated 
images of the surface brightness distribution to a ra- 
dius of 2500 AU. The computational time for a total 
5000 x 5000 x 5000 AU data cube was long, so we only 
calculated to a depth of 250 AU, keeping the field of view 
(FOV) 5000 x 5000 AU. The fluxes of the generated im- 
ages were normalized (to ensure that the geometry was 
the main constraint of the fit and not surface brightness 
variations) and rotated to angles ip = ±20° with 1° steps. 
After rotation, both the model images and the observed 
image were masked with zeros where there was no de- 
tectable surface brightness in the observed image. 

The x 2 °f the deviations of the model from the ob- 
served image were calculated. We were able to constrain 
the rotation angle of the model and the relative veloc- 
ity of the cloud to the star. The \ 2 values in the ip vs. 
i> re i phase space are shown in Figure [7] (left panel). The 
small values at large rotation angles are artifacts due to 
the masking. The best-fit rotation angle is at if = —4° to 
our initial guess, which means that the direction of mo- 
tion is 143° (CCW) of N. This is just 21° from the proper 
motion direction. The ISM velocity predicted from vec- 
torial velocity summation to fit this angle is 24 km s , 
which is close t o the ISM velocity value calculated by 
lLallement et al.l ([1995). The tangential velocity direc- 
tion of the ISM from the summation is ~ 47° CW of N, 



8 



Gaspar et al. 




Fig. 11. — Top row: Left panel: Observed 24 (im image. Center panel: Model 24 fim image including both stellar photosphere and 
bow shock. Right panel: Model image subtracted from observed. Bottom row: Left panel: Observed 70 fim image. Center panel: 
Model 70 fim image including both stellar photosphere and bow shock. Right panel: Model image subtracted from observed. The FOV is 
~ 2f7 X 2/1, N is up and E to the left. 



needed to give the observed brightness using a r ext radius 
sphere. Since the line-of-sight distribution of the dust is 
not cut off at r cxt , we used pdust = 6.43 x 10~ 26 g cm" 3 
as an initial guess; a range of density values was explored 
with model images. 

Since the surface brightness scales with the density, 
only one image had to be computed, which could be 
scaled afterwards with a constant factor. The resultant 
X 2 distribution is shown in Figure [8] The calculated 
best fitting ISM density is 5.8 ± 0.4 x 10~ 24 g cm" 3 , as- 
suming the average 1:100 dust to gas mass ratio. The 
error was calculated at A\ 2 = 1- This density (n ~ 
3.5 atoms cm -3 ) is only moderately higher than the av- 
erage galactic ISM density (~ 1 atom cm" 3 ). The calcu- 
lated surface brightness images for the three MIPS wave- 
lengths are shown in Figure [9] The closest stagnation 
point is for the 0.005 /im grains at 64 AU, while the fur- 
thest is at 227 AU for 0.056 /im grains. The temperature 
coded image in Figure \W\ shows the surface brightness 
temperature of the bow shock (i.e. the temperature of a 
black body, that would give the same surface brightness 
in the MIPS wavelengths as observed). Table CQ shows 
good agreement between the model and the measured 
values. 

The original observed images at 24 /im and 70 /im 
were compared to the model. We generated model im- 
ages with high resolution that included the bow shock 
and the central star with its photospheric brightness 
value at the central pixel. We convolved these images 
with a 1.8 native pixel boxca r smoothed STinyTim PSF 
(see lEngelbracht et all l2007t ) . These images were sub- 
tracted from the observed ones (Figure [TTj) . The resid- 
uals are small and generally consistent with the ex- 
pected noise. Finally, the best fitting SED of the sys- 
tem (r cx t = 1366 AU) is plotted in Figure [121 The total 
mass of the dust inside the r ex t = 1366 AU radius is 




1 1 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

Pfittet/Pcalculated 

Fig. 8. — The final ISM density was determined from the best 
fitting surface brightness image. This plot shows the x 2 of the fits 
of the model to the observed image, where p ca i cu iated is the initial 
guess from Figure [5] right panel and pitted i s fitted density using 
surface brightness values. 



which is pointing only 4?5 south from the galactic plane. 

The v rc \ parameter and its error are calculated by fit- 
ting a Gaussian to the phase space values at (p = —4° 
(Figure right panel) . One a errors are given by val- 
ues at Ax 2 = 1. The fits give v Tci = 35.8 ± 4.0 km s" 1 . 
Figure [5] shows that if u rc i is constrained, then we can 
also determine the density of the cloud from simple SED 
modeling. The vertical cut of Figure [5] (left panel) at 
v TC \ = 36 km s~ is shown in the right panel of the 
same Figure, pdust = 6.43 x 10" 26 g cm" 3 is derived 
from this fit, which gives an original ISM density of 
6.43 x 10 -24 g cm -3 assuming the usual 1:100 dust to gas 
mass ratio. This pdust is an upper estimate of the actual 
value, since the model computes what density would be 
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Fig. 9. — Top panels: Calculated high resolution surface bright- 
nesses for 24, 70 and 160 /im, respectively. Bottom panels: The 
24 and 70 fim image with MIPS resolution, convolved with STiny- 
Tim PSFs. The images are not rotated to the same angle as the 
observed bow shock. 
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Fig. 10. — Image (online version in color) of the bow shock gener- 
ated by the model computations. The image's FOV is 2f41 X 2'41. 
The colorscale shows the integrated surface brightness temperature 
of the bow shock (and not the radial temperature gradient of the 
grains) in Kelvins. 



Mdust = 1-706 x 10 24 g (0.023 Af M oo„). 

4. DISCUSSION 
4.1. Bow Shock Model Results 

Our model gives a consistent explanation of the total 
infrared excess and the surface brightness distribution of 
the bow shock structure at S Velorum. The question still 
remains how common this phenomenon is among the pre- 
viously identified infrared-excess stars. Is it possible that 
many of the infrared excesses found around early-type 
stars result from the emission of the ambient ISM cloud? 
The majority of infrared excess stars are distant and can- 
not be resolved, so we cannot answer for sure. However 
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Fig. 12. — The best fit SED. The window in the upper-right 
corner is a magnified part of the SED between 20 and 80 /^m. The 
plotted fluxes are 24 and 70 /im MIPS and 25 and 60 /im IRAS, 
with errorbars. The 9.7 fim silicate feature in the model SED of the 
ISM cloud is very faint and on a bright continuum. The flux from 
the G dwarf has been subtracted from the 24 and 70 fim MIPS 
observations. 



the excess at S Velorum is relatively warm between 24 
and 70 ^m (F(24) ~ 0.17 Jy, F(70) ~ 0.14 Jy), and 
such behavior may provide an indication of ISM emis- 
sion. This possibility will be analyzed in a forthcoming 
paper. Another test would be to search for ISM spectral 
features. The ISM 9.7 /im silicate feature of the dust 
grains would have a total flux of ~ 1 mJy for 6 Velo- 
rum. Since the ~ 1 mJy flux would originate from an 
extended region and not a point source that could fit in 
the slit of IRS, it would be nearly impossible to detect 
with Spitzer. Only a faint hint of the excess is visible in 
the 8 /im IRAC images, consistent with the small output 
predicted by our model. 

4.2. ISM Interactions 

To produce a bow shock feature as seen around 5 Velo- 
rum, the star needs to be luminous, have a rather large 
relative velocity with respect to the interacting ISM, and 
be passing through an ISM cloud. A relative velocity 
of ~ 36 km s~ is not necessarily uncommon, since the 
ISM in the solar neighborhood has a space velocity of 
~ 26 km s _1 ()Lallement et al.|[l9 95) and stars typically 
move with similar speeds. If the ISM encountering the 
star is not dense enough the resulting excess will be too 
faint to be detected. The Sun and its close (~ 100 pc) 
surrounding are sitting in the Local Bubble (n(HI) < 
0.24 cm" 3 , T w 7500 K. lLallementlll998t Uenkinsll2002f) . 
This cavity generally lacks cold and neutral gas up to 
~ 100 pc. The density we calculated at 5 Velorum is ~ 15 
times higher than the average value inside the Local Bub- 
ble. Observations over the past thirty years have shown 
that this void is not completely defici ent of material, but 
conta i ns filaments and cold clouds (IWennmacher et aLl 



1992t iHer bstmcicr & Wennmacher 19981 Uenkind 1200 



Mever et al.l |2006[ ). iTalbot fc Newman (|1977f) calcu- 



lated that an average galactic disk star of solar age has 
probably passed through about 135 clouds of n(HI) > 
10 2 cm -3 and about 16 clouds with n(HI) > 10 3 cm" 3 . 
Thus the scenario that we propose for S Velorum is plau- 
sible. 
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4.3. Implications for Diffusion/ Accretion Model of A 
Bootis Phenomenon 

Holwegc r et al.l ((1999) list S Velorum as a simple A 
star, not a A Bootis one. We downloaded spectra of the 
star from th e Appalachian St ate University Nstars Spec- 
tra Project (|Grav et al.ll2006l ). The spectra of S Velorum, 
A Bootis (prototype of its group) and Vega (an MK AO 
standard) are plotted in Figure [13] The metallic lines 
are generally strong for 5 Velorum. One of the most dis- 
tinctive characteristics of A Bootis stars is the a bsence 
or ex treme weakness of the Mgll lines at 4481 A (Gray 
1988). Although the Mgll line seems to be weaker than 
expected for an AO spectral type star, it still shows high 
abundance, which confirms that S Velorum is not a A 
Bootis type star (Christopher J. Corbally, private com- 
munication). The overall mctallicity ratio for 6 Velorum 
is [M/H] = -0.33 while for A Bootis it is [M/H] = -1.86 
([Gray et al.ll2006t ). The G star's contribution to the total 
abundance in the spectrum is negligible, because of its 
relative faintness. We used spectra from the NStars web 
site to synthesize a A1V/A5V binary composite spec- 
trum and found only minor differences from the AlV 
spectrum alone. Thus, the assigned metallicity should 
be valid. 




3800 3900 4000 4100 4200 4300 4400 4500 4600 



Wavelength [A] 

Fig. 13. — The spectra of <5 Velorum (bottom line), Vega (middle 
line) and A Bootis (top line). 



These results show that at 5 Velorum, where we do 
see the ISM interacting with a star, there is no sign 
of th e A Bootis phenomenon or jus t a very mild ef- 
fect. iTurcotte fc Charbonneaul (|1993h modeled that an 
accretion rate of ~ 10~ 14 M Q yr -1 is necessary for a 
T c s — 8000 K main sequence star to show the spectro- 
scopic characteristics of the phenomenon. The A Bootis 
abundance pattern starts to show at 10~ 15 Mq yr -1 and 
ceases at 10~ 12 Mq yr -1 . To reach an ISM accretion of 
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10 -15 Mq yr -1 a collecting area of 2 AU radius would 
be needed with our modeled ISM density and velocity. 
For an accretion of 10~ 14 , 10~ 13 and 10~ 12 M yr" 1 
collecting areas of 6.5, 20 and 65 AU radii are needed, 
respectively. 

With the accretion theory of lBondi fc Ho vlc (1944), we 
get an accretion rate of 6.15 x 10~ 15 Mq yr -1 for 5 Velo- 
rum. Thus, the accretion rate for this star is probably 
not high enough to show a perfect A Bootis spectrum, 
but should be high enough for it to show some effects of 
accretion. This star is an exciting testbed for the diffu- 
sion/accretion model of the A Bootis phenomenon. 

5. SUMMARY 

We observe a bow shock generated by photon pressure 
as 5 Velorum moves through an interstellar cloud. Al- 
though this star was thought to have a debris disk, its 
infrared excess appears to arise at least in large part from 
this bow shock. We present a physical model to explain 
the bow shock. Our calculations reproduce the observed 
surface brightness of the object and give the physical pa- 
rameters of the cloud. We determined the density of the 
surrounding ISM to be 5.8 ± 0.4 x 10 -24 g cm" 3 . This 
corresponds to a number density of n s=s 3.5 atoms cm~ 3 , 
which means a ^ 15 times overdensity relative to the av- 
erage Local Bubble value. The cloud and the star have 
a relative velocity of 35.8 ± 4.0 km s~ . The velocity of 
the ISM in the vicinity of 5 Velorum we derived is consis- 
tent w ith LIC velocity measurements bv lLallement et all 
(1995). Our best-fit parameters and measured fluxes are 
s ummarized in Ta bic [1] 

iHolweger et al.l ()1999[ ) found that d Velorum is not a A 
Bootis star. The measurements from the Nstars Spectra 
Project also confirm this. Details regarding the diffu- 
sion/accretion time scales for a complex stellar system 
remain to be elaborated. Nevertheless, our Spitzer ob- 
servations of 8 Velorum provide an interesting testbed 
and challenge to the ISM diffusion/accretion theory for 
the A Bootis phenomenon. 
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